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c^ ■ Abstract 

Dirac monopoles have been widely studied and searched, though never found. A 
way out of this impasse is the idea that monopoles are not seen freely because they are 
confined by their strong magnetic forces forming a monopole-antimonopole bound state 
called monopolium. Monopolium was shown to be produced abundantly and in some 
scenarios easier to detect than monopoles themselves. The Large Hadron Collider is 
reaching energies never achieved before allowing the search for exotic particles in the TeV 
mass range. We extend our previous analysis to the observability of monopolium at LHC 
in the 7 7 channel particularizing our quantitative discussion to monopolium masses that 
can be detected with integrated luminosites ~ 1 fb~^. 
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1 Introduction 

Numerous experimental searches for magnetic monopoles have been carried out but all 
have met with failure p1-[T0]. These experiments have led to a lower mass limit in the 
range from 350 — 500 GeV. 

This lack of experimental confirmation has led many physicists to abandon the hope 
in their existence. A way out of this impasse is the old idea of Dirac [TTHT3] . namely, 
monopoles are not seen freely because they are confined by their strong magnetic forces 
forming a bound state called monopolium [T31[T2]. This idea was the leitmotiv behind our 
research, namely we proposed that monopolium, due to its bound state structure, might 
be easier to detect than free monopoles p^[T7] . We showed that certain parameterizations 
of the mass and the width, allowed for such scenario. 

The Large Hadron Collider (LHC), which entered last year in operation colliding 3.5 
TeV protons, will probe the energy frontier opening possibilities for new physics includ- 
ing the discovery of magnetic monopoles either directly, a possibility contemplated long 
time ago [18l[T9], and revisited frequently [3|[7t [201425] . or through the discovery of mo- 
nopolium, as advocated in refs. [T6l[T7]. The direct observation is based on the search 
for highly ionizing massive particles at the ATLAS [26] and CMS [27] detectors or at 
the MoEDAL experiment ^28j|, the design of which is optimized to search precisely for 
such exotic states. These developments motivate our present research which analyzes the 
production of monopolium at LHC by the mechanism of photon fusion and its subsequent 
decay into 77. 

As was shown in our previous work, photon fusion is ideal for monopolium production 
[T7] . To complete the analysis, we study here how monopolium, produced at the collision 
point, can be detected through its decay in a pair of extremely energetic photons, a 
channel for which LHC detectors are optimized, since this is one of the decay channels 
for finding a light Higgs with low counting rates. 

In the next section we establish the formalism for monopolium production and decay. 
Section 3 describes how to incorporate the elementary process into p—p scattering. Section 
4 presents our results and studies the variation with the two parameters, monopole mass 
and monopolium mass. Finally in section 5 we draw some general conclusions of our 
study. 

2 Monopolium dynamics 

For the case of monopole interactions at energies higher than their mass there is no 
universally accepted effective field theory [2UH^ . However, the study of the classical 
interaction of a monopole passing by a fixed electron leads to an interaction for the 
monopole which is associated with the electric field felt. If one compares this interaction 
with that of a positron passing by an electron, one realizes that the difference between 
QED and dual QED [31], the theory without strings for monopole interactions, is simply 
to change the electric charge by the magnetic charge times the velocity. For this reason we 
will employ a minimal model of monopole interaction which assumes an effective monopole 



photon-coupling which is proportional to the nionopole's induced electric field gP for a 
monopole moving with velocity /3 [5| [2T1 - [23] . 

The Dirac quantization condition does not specify the spin of the monopoles. We 
choose here monopoles of spin 1/2 coupled in monopolium of spin in order to have a 
minimum energy radial structure. 




Figure 1: Diagrammatic description of the elementary subprocess of the monopolium production 
from photon fusion. V{r) represents the interaction binding the monopole-antimonopole pair to 
form monopolium. 

Recently we studied the production of monopolium by photon fusion at LHC [TTj. 
The elementary subprocess calculated is shown in Fig. [1] The standard expression for 
the cross section of the elementary subprocess for producing a monopolium of mass M is 
given by 

where T{E), with E off mass shell, describes the production cross section. Note that 
r[M] = 0. T M arises from the softening of the delta function, 5{E'^ — M^) and therefore 
is, in principle, independent of the production rate T[E) and can be attributed to the 
beam width |32l|33]. 

In Fig. [2] we show the total cross section for monopolium production from photon 
fusion under present LHC running conditions for a monopole mass (m) ranging from 500 
to 1000 GeV. In the figure the binding energy is fixed for each mass (2 ?ti/15), chosen so 
that for our case study, m = 750 GeV, the binding energy is 100 GeV and thus M = 1400 
GeV. With this choice the monopolium mass (M) ranges from 933 to 1866 GeV. We 
notice that detection would be possible with an integrated luminosity of 1 fb~^ if the 
chosen binding energy is at the level of 10% of the monopole mass or higher. In the 
present analisis we study binding energies small compared to the bound state mass, M, 
in order to be consistent with the formalism used. 
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Figure 2: Total cross section for monopolium production at LHC with 3.5 TeV beams for 
monopole masses ranging from 500 to 1000 GeV (full curve). The broken curves represent the 
different contributions to the total cross section as described in section 3: semielastic (dashed), 
elastic (shorst dashed) and inelastic (dotted). We have chosen a binding energy ~ 2 m/15 and 
Vm = 10 GeV. 



The interest in this paper is in the detection of photons after monopolium decay. The 
elementary subprocess is shown in Fig. [3l which could be considered as a contribution to 
light-by-light scattering in the presence of monopolium. 

The standard expression for the cross section of this elementary subprocess, after 
having integrated over angles, is given by 



cr(77 — )■ M — )• 77) 
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E^ (E^-MY + M2 Tlj 
Here T{E), with E off mass shell, describes the vertex 77M. Monopolium is stable in 
the center of mass but we add an experimental Gaussian width Tm ~ 10 GeV in line with 
the values used in ref. [341. 
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Figure 3: Diagrammatic description of the monopolium production and decay. 



We recall now the computation of the T{E), which represents the vertex of the mo- 



nopolium decay to 7 7 . The calculation, following standard field-theory techniques of the 
decay of a non-relativistic bound state, leads to 
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We have used the conventional approximations for this calculation: the monopole and 
antimonopole, forming the bound state, are treated as on-shell particles, when calculating 
the elementary scattering process shown on the right of Fig. HI the bound state is described 
by a wave function obtained from a Coulomb- type interaction between the pair [T7 tl32ll33] . 
However, once the calculation is performed we substitute 2m by M, where m is the 
monopole mass to take into account binding. In the expression, ag corresponds to the 
photon-monopole coupling and ipM is the monopolium ground state wave function. 




Figure 4: The monopolium vertex and its microscopic description . 



Before we proceed let us raise a warning on our calculation. Bound relativistic 
systems are notoriously difficult dynamical objects. We proceed here by performing a 
non-relativistic calculation of the monopolium wave function. The validity of the non- 
relativistic approximation in the bound state wave functions was analyzed in a previous 
work [16]. In the dynamics of the decay formula the substitution of 2m by M constitutes 
an intuitive way to take the off-shellness of the monopoles, i.e. their binding energy in 
monopolium, into account. For the purposes of estimation both approximations seem 
reasonable specially since our binding energies will never exceed 15% of the monopole 
mass, i.e. less than 10% of the total mass of the system. 

Using the Coulomb wave functions of ref. [16] expressed in the most convenient way 
to avoid details of the interaction, which will be parameterized by the binding energy, one 
has 
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and the effective monopole coupling theory described above in the case of monopolium 
production, gives rise to [T7] . 
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Figure 5: Elementary cross section as a function of w = E/M in units of 1/M^ calculated for 
m = 750 GeV and M = 1400 GeV. Left: near threshold. Right: away from threshold. 
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Here, a is the fine structure constant and /3 the monopohum velocity, 



/5 




(6) 



which is the velocity of the monopoles moving in the monopolium system. 

Note that due to the value of /3 the vertex vanishes at the monopolium mass, where 
the velocity is zero. Therefore a static monopolium is stable under this interaction. We 
refer to refs. [5l[23] for a thorough discussion on Lorentz invariance of the theory. 

A caveat is due here. There is a duality of treatments in the above formulation as can 
be seen in Fig. [H The static coupling is treated as a Coloumb like interaction of coupling 
g binding the monopoles into monopolium, although ultimately the details are eliminated 
in favor of the binding energy parameterized by the monopolium mass M. We find in this 
way a simple parametric description of the bound state. The dynamics of the production 
of the virtual monopoles, to be bound in monopolium, is described in accordance with 
the effective theory [2n[23] . and this coupling is Pg. This is similar to what is done 
in heavy quark physics [22] (see his figure 5), where the wave function is obtained by a 
parametric description using approximate strong dynamics while the coupling to photons 
is elementary. 

The production cross section can now be written as. 
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The above cross section satisfies comfortably the unitarity limit [36] , 

TT 



a < 



3^2- 



p{X) 




Figure 6: Processes contributing to the 77 cross section. The blob contains the three 
cases described in the text, p represents the proton and X any other hadronic state. 



To feel safe with our approximations we consider the binding energy much smaller 
than m, i.e. M ~ 2m. In this case the elementary cross section has two very different 
behaviors as shown in Fig. [5] : i) at threshold it is dominated by /3 and the cross section 
rises (see left figure); ii) away from theshold the dominant behavior is the 1/E dependence 
and the cross section drops faster than the unitary limit. The conflict between these two 
behaviors produces a wide bump-like structure. 

3 Analysis of p-p scattering 

LHC is a proton-proton collider, therefore, in order to describe the production and desin- 
tegration of the monopole-antimonopole pair, we have to study the following processes 
above the monopole threshold {(3 > ^^ E > M), 



p + p -^ p(X) +p(X) + 7 + 7, 



(9) 



shown globally in FigJHl where p represents the proton, X an unknown flnal state and we 
will assume that the blob is due exclusively to monopolium. This diagram summarizes 
the three possible processes: 

i) inelastic p + p->X + X + 7 + 7->X + X + M^X + X + M + 7 + 7 
ii) semi-elastic p + p^p + X + ^y + 'y^p + X + M-^p + X + M + 'j + 'j 
in) elastic p + p-^p + p + y + y^p + p + M^p + p + M + 'j + 'y. 

In the inelastic scattering, both intermediate photons are radiated off partons (quarks 
or antiquarks) in the colliding protons. 



In the semi-elastic scattering one intermediate photon is radiated by a quark (or 
antiquark), as in the inelastic process, while the second photon is radiated from the 
other proton, coupling to the total proton charge and leaving a final state proton intact. 

In the elastic scattering both intermediate photons are radiated from the interacting 
protons leaving both protons intact in the final state. 

In the blob we incorporate the elementary subprocess shown in Fig. [3] and described 
by Eq. (JH). 

We calculate the 77 fusion for monopolium production following the formalism of 
Drees et al. |37l|38]. 

In the inelastic scattering, p + p— ;'X + X + 7 + 7— )'X + X + M + 7 + 7, to approx- 
imate the quark distribution within the proton we use the Cteq6-1L parton distribution 
functions [39] and choose Q^ = s/4 throughout, where s is the center of mass energy of 
the elementary process. 

We employ an equivalent-photon approximation for the photon spectrum of the in- 
termediate quarks jlHlllI]. 

In semi-elastic scattering, p + p— >p + X + 7 + 7— )-p + X + M + 7 + 7, the photon 
spectrum associated with the interacting proton must be altered from the equivalent- 
photon approximation for quarks to account for the proton structure. To accommodate 
the proton structure we use the modified equivalent-photon approximation of [38] . 
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Figure 7: The differential cross section (solid line) and its components (semielastic 
(dashed), elastic (short dashed) and inelastic (dotted)) as a function of the center-of- 
mass photon energy for a monopole mass of 750 GeV and a monopolium mass of 1400 
GeV. 



The total cross section is obtained as a sum of the three processes. The explicit 
expressions for the different contributions can be found in [22] . 

In order to obtain the differential photon-photon cross section from the above formal- 
ism we develop a procedure which we exemplify with the elastic scattering case. In that 
case the p — p cross section is given by [38] . 



o-pp{s) = dzi I dz2f{zi)f{z2)(ryy{ziZ2s), (10) 

JM'^/s JM'^/szi 

where -y/i is the center of mass energy of the pp system and the /'s represent the elastic 
photon spectrum. 

We perform the following change of variables 

V = Z1Z2 , w = Z2 , 

which leads to 



app[s) 



dv —f{-)f{nj)cr,,{vs). (11) 

JaP/s Jv w w 

Note that to fix the center of mass energy of the photons is equivalent to fix v. For fixed 
V we have, 

which can be rewritten in terms of E^, the center of mass energy of the photons, and the 
elementary photon-photon cross section as, 

^(E,) = ^ a,,{s,,) f —f{^)f{w). (13) 

This procedure can be generalized easily to the semielastic and inelastic cases, where the 
appropriate change of variables are 

V = Z1Z2X1 , W = Z2X1 , U = Xi 

and 

V = Z1Z2X1X2 , W = Z2X1X2 , U = X1X2 , t = X2, 

respectively which one has to introduce into integral expressions with a product of three 
/'s (semielastic) or four /'s (inelastic) representing quark densities and photon spectrum 

[371EH]. 

4 Results and Discussion 

Our aim is to show scenarios which could arise during the present LHC running period 
and to discuss general properties of the monopolium system which might serve when 
higher luminosities are achieved. In Fig. [7] we show the structure for the differential cross 
section. It is a wide bump, starting very close after threshold, i.e. the monopolium mass 
(1400 GeV in this case), and extending for about 1000 GeV . We show in the figure the 
contribution of the different components to the cross section. The elastic and semielastic 
components dominate. The behavior is well understood by the structure of Eq. [TJ the 
bump initiates due to the rising of the cross section close to threshold associated with 
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Figure 8: Differential two photon cross section as a function of tlie pfioton center-of-mass 
energy for a monopole mass of 750 GeV and different binding energies for monopolium: 
solid line 75 GeV, dasfied 150 GeV and solid 225 GeV. 



its (3 beliavior. Close to threshold /3 takes almost its asymptotic value of 1 and the 1/E 
behavior of the cross section starts to softly dominate (recall Fig|5]). 

Two are the main physical dependences of the cross section: the monopole mass, m, 
and the monopolium mass , M. In Fig. |8] we fix the monopole mass to 750 GeV and 
vary the binding energy. We see that the cross section increases dramatically with binding 
energy. Thus the effect of the binding is twofold: it increases the cross section and it lowers 
the threshold from 2m the monopole-antimonopole production threshold. To observe the 
other dependence, in Fig. |9]we fix the binding energy to 100 GeV and vary the monopole 
mass. The effect goes inversely proportional to the monopole mass, i.e., the lower the 
monopole mass the larger the cross sections. 

We can summarize our findings by stating that low monopole masses and large bindings 
favor the detection of monopolium. Monopolium has the advantage over the monopole- 
antimonopole process studied in ref [21] of lowering the threshold, narrowing the bump 
and increasing the cross section with binding energy. In Fig. [10] we show three interesting 
effects, namely the two photon decays of Higgs, monopolium and monopole-antimonopole 
annihilation. The Higgs signal has been increased over the background by a factor of 
50. The parameters for the monopolium cross section shown are m = 750 GeV and 
M = 1400 GeV, which have been chosen so that its signal is of the same size of that 
of the monopole-antimonopole annihilation as calculated in ref. [21] with a monopole 
mass of also 750 GeV. We note two of the features mentioned before, the lower threshold 
and the narrower bump structure. If we would increase the binding by a few tens of 
GeV the height of the bump would increase considerably with respect to the monopole- 
antimonopole cross section (recall FigJHj). Note that the 77 conventional background at 
the monopole scenario is extremely small, as measured recently by ATLAS in di-photon 
studies [12] and while searching for ilf — )■ 77 |43] . 
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Figure 9: Differential cross section as a function of photon energy for different monopole 
masses: 500 GeV (solid), 750 GeV (daslied), 1000 GeV (dotted) and fixed binding energy 
(100 GeV). 



If a broad bump would appear experimentally one should confirm that it could arise 
from monopole dynamics. Thereafter the way to distinguish between the decay of mo- 
nopolium and the annihilation of monopole-antimonopole would be through the angular 
dependence. In the present scenario, a spin zero monopolium, the angular decay properties 
would be similar to that of para-positronium [15], while that of monopole-antimonopole 
would be analogous to light-by-light scattering in QED [l^HT]. Moreover, the two pe- 
henomena could occur simultaneausly, as happens in the case of electrons and positrons, 
where we have in light-by-light scattering electron-positron annihilation and positronium 
decay. If the latter were the case two bumps, if the overlap is not large, or a very broad 
fiat bump, if the overlap is considerable, could be seen. The existence of one or two bumps 
depends very strongly on the binding dynamics and the monopole mass. Note that there 
is no possible confusion with the Higgs, since its width is narrow compared to its mass. 
Moreover, the scattering cross section for a heavy Higgs in the two photon channel is 
extremely small compared to that in the other channels, and therefore, its characteristics 
would be known by the time two photon bump would be seen. 



5 Conclusions 

The Dirac quantization condition is a beautiful consequence of the existence of monopoles 
and therefore they represent an extremely appealing physical scenario. There is as of yet 
no experimental proof of their existence. This has led to approximate mass bounds which 
suggest a mass scale for the monopole above 500 GeV. LHC opens up this energy regime 
for research and therefore monopoles become again a subject of experimental search. 
Even if monopoles exist it might be possible that due to the very strong magnetic 
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Figure 10: We show the differential two photon crossection with the Higgs signal (scaled 
by a factor of 50 over the background) where the monopolium and monopole-antimonopole 
contributions have been incorporated. We also show the monopolium signal (dotted) and 
monopole- antimonopole as in ref. [21] (dashed) (scaled by a factor of 5). The mass of 
the monopole in both cases is 750 GeV and that of the monopolium 1400 GeV. Note that 
for Ifb^^ the y-axis measures counts in 1 GeV bins. The background has been obtained 
from ref. 




coupling they do not appear as free states but bound forming monopolium, a neutral state, 
very difficult to detect directly. Here we have analyzed the coupling of monopolium to 
photons and its contribution to light-by- light like scattering, i.e. monopolium's dynamical 
decay. We have found that for reasonable values of the monopole mass and relatively small, 
compared with their mass, binding energies, spin zero monopolium desintegrates into two 
gammas with cross sections which are reachable in this first run of LHC. 

Our investigations go beyond this wishful scenario. We have seen that the cross 
section depends both on the mass on the monopole and on that of monopolium, increasing 
inversely proportional to the monopole mass and being directly proportional to the binding 
energy. This means that similar cross section can be achieved with very heavy monopoles 
if the binding energy is large. This is however a qualitative statement since our calculation 
has been carried out in a non relativistic framework for which large bindings are not well 
under control. 

To conclude, the eventual appearance of broad bumps in the two photon cross section 
might be associated with the existence of monopoles, either free, as discussed in our 
previous work [21j, and/or bound in monopolium as presented here. The implementation 
of our findings in detector analysis would provide actual observations. 
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